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Abstract 

We have studied the distribution of equihbrium edge current density in 2D system in 
a strong (quantizing) magnetic field. The case of half plane in normal magnetic field has 
been considered. The transition from classical strong magnetic field to ultraquantum 
limit has been investigated. We have shown that the edge current density oscillates and 
decays with distance from the edge. The oscillations have been attributed to the Fermi 
wavelength of electrons. The additional component of the current smoothly depending 
on the distance but sensitive to the occupation of Landau levels has been found. The 
temperature suppression of oscillations has been studied. 

The magnetic field acting on the low dimensional system is traditionally considered as 
homogeneous and coinciding with the external field. Nevertheless, the magnetization of 
such quantum system as atom substantially changes the magnetic field in the vicinity of 
nucleus. This is essential for some phenomena, e.g., NMR. The distribution of magnetic 
field in the systems with spatial quantization have been studied by authors in papers [T] 
using the linear approximation in an external field. This problem is close to the problem 
of orbital magnetism intensively investigated in different systems with separable and non- 
separable variables |21 El E! • It was shown, that the susceptibility of a large system at low 
temperature strongly fluctuates and changes the sign when the Fermi level moves between 
the energy levels of the system. 

The goal of the present paper is to study the spatial distribution of equilibrium currents 
in 2D semi-infinite plane sample subjected to a strong magnetic field. Such mathematical 
setting is applicable for a sample of arbitrary shape if its characteristic dimensions exceed 
the inverse cyclotron diameter. 

Edge current in 2D system in a finite magnetic field. 

Let us consider a semi-infinite plane x>0, — oo < y < oo, with hard border in a magnetic 
field Bz = B. Assume the vector potential has gauge Ay = Bx. The states in the presence 
of magnetic field can be described by the longitudinal momentum p and the transversal 
number n: ipn, p{x)e^^y. Here and below h = 1. The wave functions can be expressed via the 
parabolic cylinder function D^{x): 

oo 

V'n,p(x) = CD,^ [V2 {x - Xp) /a) , J D,^ [V2x/a) dx = 1. (1) 
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Here a = yc/eB is the magnetic length, xq = —pa^. The boundary condition Di,^{—\/2pa) = 
determines the energy levels En,p = co'(z/„(p) + l/2), n = 0, 1...; u = eB/nieC is the cyclotron 
frequency, rUe is the electron effective mass. 
The density of current has the form 



where jo = eu/a = e/{mea^) is the characteristic current density, f± = f{En{p) ± gfisB /2), 
f{E) = (exp {{E — n)/T) + 1)^^ is the Fermi distribution function (/i, T are the chemical 
potential and temperature, correspondingly), g is the (^-factor. The expression (2) contains 
two contributions caused by the orbital and spin parts of the current (the first and the second 
lines, correspondingly). Below we shall neglect the spin term, assuming the smallness of the 
(7-factor. 

Plots of the edge current calculated according to Eqs.(IT)) and Q are presented in Figs. 1-3. 
Let us discuss behavior of the edge current at low temperature. The direction of the current 
is determined by the vector product of the normal and the magnetic field. However, it does 
not mean the constancy of the sign of the surface current density. Really, let only one Landau 
level be occupied. Consider the states localized far from the boundary. These states are not 
perturbed by the wall. The contribution to the current density from one state with given p 
is antisymmetric relative to the point Xp. As the state approaches to the boundary the level 
goes up. Occupied levels (they lie under Fermi level) have close wave functions. Averaging 
over momentum compensates the current density far from the point of intersection the Fermi 
level with Landau level, while near this point noncompensated contribution of constant sign 
(positive in this instance) remains. When the Fermi level approaches the lowest Landau 
level, the distance of the edge state from the border grows as ln^/2(/i/cj - 1/2) and the edge 
state becomes more and more ideal. Under increase of Fermi level the wave function begins 
to distort due to boundary, positive part of the contribution to the average current density is 
suppressed but the negative contribution appears and the graph of current density acquires 
negative minimum. 

Crossing next Landau levels results in new contributions to the current density. First 
these contributions are situated on large distance from the boundary, then they move to it 
and merge with contributions from lower lying states. The upper is the state the wider is the 
region it occupies and the more oscillations it contributes to the current density. When we 
have many occupied Landau levels, the contributions from different levels merge to produce 
the Friedel oscillations of surface current density studied by authors earlier [T] . However, the 
number of oscillations is restricted (unlike earlier considered limit of weak magnetic field), 
because the surface current is distributed on the thickness of the order of cyclotron diameter 
2rc. Within this region the current density oscillates, on the outside it decays exponentially. 
The number of oscillations is defined by the number of occupied Landau levels. In the region 
X <^ IVf. the effect of magnetic field is weak and the linear-response expression from pQ is 
valid. 

As it is seen from Fig. 2 the smooth dependence on coordinate x is superimposed on 
the space oscillations. The total edge current J experiences the alternating-sign oscillations 
with magnetic field. These oscillations accord with the de-Haas-van Alphen oscillations of 
magnetic moment. 



jo ( 

jyi^) = J dpali^p - ^n,p(a;)(/+ + /-) 




(2) 
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Note that in the hmit of weak magnetic field the smooth contribution transforms to 
studied in |lj contribution depending hnearly on the distance, and oscillations decaying on 
power law spread to infinite range from the boundary. 

Temperature suppresses the oscillations of current density (see Fig. 2). It happens more 
effectively on large distance from boundary. The suppression takes place due to temperature 
dephasing of electrons near the Fermi surface on characteristic length 1^ = kp/{27TmeT), 
where kp is the Fermi momentum pp. At finite magnetic field the competition of two lengths 
which cut down oscillations (cyclotron diameter and It) occurs. The scattering suppress both 
spatial oscillations and smooth contribution as well. 

The alternating-sign oscillations of total edge current seem at first sight strange if to take 
into account, that at large Fermi energy the total edge current should be diamagnetic and 
be described by the formula 



12n 

However, in the limit of the large system the total quantity of edge current J is directly 
connected to the magnetic moment of the system by the relation M = JS/ c, where S is the 
system area. At the same time, the total moment at T = can be found using f2-potential 
of the system 



M = -— = ^^^^ V 
dB nc to 



--(2n + l) 



e[li-[n + ^]u]. (3) 



This expression oscillates with /i, jumping when the Fermi energy crosses the Landau levels. 
In accord with this formula the moment vanishes when the Fermi energy lies in the middle 
between the Landau levels. 

On the other hand, the edge currents in nonequilibrium conditions are connected with 
the quantization of microcontact resistance in adiabatic transport regime of quantum Hall 
effect. The total current in the state {n,p) is given by expression 

3n, = -e^. (4) 

As the states are localized on x, the edge current is determined by the states with Xp, lying 
near to the border. The partial current density j„ p exponentially decreases with increase of 
Xp. Therefore, summing currents over all states near the border we obtain finite quantity. 
Such reasons determine the quantization of microcontact resistance. 
The total edge current in semi-infinite problem at T = looks like 



'/ = -eE%^(/^-^n(p)) = -^E 

n,p " n=0 



n -\ — 



2/ u 



9[fi~[n + ^)u;). (5) 



The current through the structure is expressed via the difference of the edge currents 
corresponding to chemical potentials of edges /xi and /i2. The difference /xi — /i2 in a nonequi- 
librium problem coincides with the potential difference V, applied to the microstructure. 
From here we conclude, that the nonequilibrium current is 

J = e^VN/h. 

The chemical potentials of edges are equal in equilibrium, but the edge current from 
one edge should be determined by the same expression, as in non-equilibrium. However, 
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such understanding of the Eq. (0) contradicts Eq. Q for the magnetic moment and numerical 
calculations. 

To find out the nature of the discrepancy, we shall find the contribution to current density 
for the boundless problem from Landau levels with a momentum lying between pi and p2, 
where \xp-^ — XpJ ^ a 

P2 . P2 pi 

jy{.Pi,P2;x) = -^^y {x-Xp)<^l (^ ^ J^^ ^ dp= -^^(J dp... - y dp...). (6) 

" Pl " -oo -oo 

Here (fniO dimensionless normalized functions of the harmonic oscillator. Contributions 
from regions close to the points Xp^ and independent and separated in the space. 

Two edge currents are presented in Eq.(jni)(and can be isolated), but cancel each other in the 
total current. The resulting contribution to a total edge current from current density near 
the point a;^^ is located far from border and should be subtracted from a total edge current. 

This part of a current is determined by integration of jy {pi ,P2',x) with respect to x near 
the point Xp 




Here A a\fN . Subtracting (|7j) from (0) we obtain the correct expression for a total edge 
current J = Mc/S, with M from the 0. Numerical calculations (see Fig. 3) accord with 
this formula. 

Conclusions. 

The presence of edge current results in weak change of a magnetic field applying to the 
2D system. The non-uniformly distributed current creates the non-uniformly distributed 
magnetic field. 

Note that the spatial inhomogeneity of a magnetic field can effect on any responses 
sharply dependent on a magnetic field, in particular, on the Shubnikov oscillations, geomet- 
rical resonances or magnetic focusing. One can expect that the inhomogeneity will result 
in washout of sharp singularities in these quantities. This gives the way to measure the 
magnetization. 

Other variants can be based on sensitivity of nuclear spins to a local magnetic field. 
NMR-markers, placed on the certain atomic planes, can act as indicators of local magnetic 
fields. One can propose to use the constant in time magnetic field for shift of NMR line 
or alternating magnetic field for excitation of transitions. The temporal change of inhomo- 
geneous magnetic field can be produced by control of the electron wave functions due to 
alternating gate voltage. As another way one can suggest the periodic change of tempera- 
ture of electron gas by source-drain voltage with subsequent change of 6B) caused by the 
temperature suppression of the Friedel oscillations. 
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Figure 1: On the left: rehef of the edge current density in units jo at finite magnetic field as 
function of the distance from the border and Fermi energy. The grey color corresponds to 
zero density of a current, white - positive, black - negative. On the right: dependence of the 
current density on distance at the selected Fermi energy (marked by straight lines in the left 
figure) . 
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Figure 2: Dependence of density of edge current on temperature. The Fermi level is close 
to the 5th Landau level: solid curves - Ep = AAuj, dotted curves - Ep = 4.6cj; in these 
cases the integrated currents are oppositely directed. At temperatures T ^ \Ep — 4.5a;| the 
curves merge. At T ^ a; the current density ceases to depend on a magnetic field. 




Figure 3: The total edge current as function of a magnetic field at T = calculated from 
Eq.(2) (dotted) and from Eq. (3) (solid). 
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